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Abstract 

Axions in the Peccei-Quinn (PQ) mechanism provide a promising solution to the strong CP 
problem in the standard model of particle physics. Coherently generated PQ scalar fields could 
dominate the energy density in the early Universe and decay into relativistic axions, which would 
conflict with the current dark radiation constraints. We study the possibility that a thermal 
inflation driven by a 17 ( 1 ) gauged Higgs field dilutes such axions. A well-motivated extra gauged 
17 ( 1 ) would be the local B — L symmetry. We also discuss the implication for the case of U(1)b~l 
and an available baryogenesis mechanism in such cosmology. 
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I. INTRODUCTION 


Inflationary cosmology has been the paradigm in modern cosmology [l]. The accelerating 
expansion in the early Universe not only solves the flatness and horizon problem jsj] but also 
provides the density perturbations which grow to the large-scale structure {sj]. The false 
vacuum energy inducing quasi-de Sitter expansion during inflation is transferred into the 
energy of radiation with high temperature through reheating processes, so that the filling of 
early Universe with high-density and high-temperature thermal plasma is realized. 

It is often expected that such a hot universe directly connects to the hot big bang universe 
in which big bang nucleosynthesis (BBN) takes place after one second. However, there are 
various possibilities of nontrivial evolution scenarios where a stage dominated by temporal 
matter or vacuum energy would exist between the primordial inflation and the onset of BBN. 
In addition, through the evolution of the early Universe from a very high-energy to low- 
energy region, it is supposed that the Universe has undergone a series of cosmological phase 
transitions and the symmetry breaking of particle theory, e.jm a large-gauge (sub)group of 
grand unified theory (GUT), Peccei-Quinn (PQ) symmetry 4|, and the electroweak gauge 
symmetry in the standard model (SM) of particle physics. There are two well-motivated 
symmetries expected to be broken at an intermediate scale from a particle physics viewpoint. 
One is a PQ symmetry introduced to solve the strong CP problem in the SM. The other is 
B — L symmetry which might be gauged symmetry at a high-energy scale and its breakdown 
would generate Majorana masses of right-handed (RH) neutrinos to account for neutrino 
oscillation phenomena Q . 

A scalar held to break the PQ symmetry has small couplings to SM particles suppressed 
by the PQ breaking scale f a , also known as the axion decay constant. The PQ scalar held 
with a “wine bottle” potential could be produced as coherent oscillation due to its scalar 
nature and temporally dominate the energy density of the Universe if its decay rate is very 
small because of suppressed couplings. The radial direction of the PQ scalar 1 would mostly 
decay into axions or SM particles through loop processes. For the KSVZ axion model 6] 
where heavy quarks carrying a PQ charge are indeed heavy enough, this could happen. For 


the DFSZ axion model 


7] where axion multiplets have direct couplings with the SM Higgs 


1 From now on, we simply call it the PQ scalar. 
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field, this could be the case if its coupling between the PQ field and the SM Higgs field is 
small enough so that the PQ radial scalar decay into Higgs bosons is suppressed. Thus, 
associated with the breakdown of the global U(1)pq symmetry, the PQ scalar likely decays 
into mostly axions in this model and cannot reheat the Universe well. Those overproduced 
massless axions act as dark radiation which is nowadays stringently constrained. The Planck 
data show the effective number of neutrino spices 8] 


N eff = 3.15 ± 0.23, 


( 1 ) 


and the difference from the SM prediction is the so-called A N e ff. 

This dark radiation axion is an unwanted relic and needs to be diluted anyway. Thermal 
inflation is a well-known mechanism to dilute such unwanted relics { 9 ] and is driven by 
a scalar held p, often called the baton, whose scalar potential V(p) is of the symmetry¬ 
breaking type with a large VEV. In a hot Universe after the reheating of the primordial 
inflation, the scalar held (p is thermally excited and stays at the origin p = 0 because the 
thermal correction in the effective potential makes the origin of the potential to a temporal 
minimum with its large false vacuum energy U(0). When the false vacuum energy becomes 
comparable to the background radiation energy density, the Universe starts accelerating 
expansion again. Since the radiation energy density is redshifted away as its temperature 
decreases exponentially, the origin tp = 0 becomes unstable and the inflation terminates. 
With the enormous entropy production by the p decay, unwanted relics can be diluted 
away. 

In this paper, we investigate the condition of successful thermal inflation to dilute axions 
generated by the late decay of the dominated PQ scalar held. If this baton p is a gauge singlet 
and has an (approximate) global (7(1), then the axions associated with the baton could be 
produced again as shown in Ref. [lOj]. Thus, in order not to have an axion overproduction 
problem after thermal inflation, we assume that a baton held is charged under a local (7(1) 
symmetry. We also discuss the implication in the case that this local (7(1) symmetry is 
identihed with gauged U{1)b~l- 
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II. THERMAL INFLATION IN AN AXION-DOMINATED UNIVERSE 


A. Thermal inflation by a wine bottle potential 


We describe the outline of the scenario with a thermal inflation period driven by canonical 
scalar held <p, which is responsible for breaking its local U( 1) symmetry. For illustrative 
purposes, in this section, we assume the scalar potential for thermal inflation is of the wine 
bottle form. The scalar potentials for tp are given by 

V(<p) = A,(u 2 - M 2 ) 2 . (2) 


A flaton held p is assumed to be in thermal equilibrium through interactions with particles 
in the hot thermal bath and, hence, the thermal correction; in fact, the thermal mass term, 


SV = ffr 2 M 2 , 


(3) 


with T being the temperature of the thermal plasma, is added in the scalar potential. Here, 
g v is parametrizing the coefficient, while sometimes we may use an effective coupling with 
another particle h = yfg^ instead of g ip in the rest of this paper. 

Before evaluating the number of e-folds in a relativistic axion-dominated Universe, we 
here review how to estimate it in the standard radiation-dominated Universe. When the 
false vacuum energy V (0) dominates over the energy density of the background radiation 
p = 7r 2 g*T 4 /30 with g* being the number of relativistic degrees of freedom, a thermal inflation 
starts. We define the initial temperature of a thermal inflation by 

= V(0). (4) 


A thermal inflation terminates when the temperature drops below the critical temperature 
T c defined by 


9 C p r7-,2 \ 2 

f 4 - T C = V > 


(5) 


and the origin p = 0 becomes unstable. In the rest of this paper, we use the final temperature 
Tf instead of Tc- We And the initial and final temperature, and the number of e-folds as 

.2 \ -i/4 


7T 


T i ( 3() . r b 


(V 4 ) 


4\ 1/4 


T f = 4i 


I3\ 

9<p 


-v, 


( 6 ) 

(7) 


4 



( 8 ) 


iV 4 = — In 4\/3 — — In 


7 

30 


9* 


- 7 In 

4 4 


Here, the subscript 4 in iV 4 indicates, for later convenience, that the scalar potential is 
quartic. We find several number of e-folds by thermal inflation can be realized only for 


-V ^ 9ipi 


from Eq. (JHJ) - One may find that g, P would be expressed as 


( 9 ) 


9<p ~ K + 9 2 + ^2 y 2 


( 10 ) 


where g and y stand for the gauge coupling of (p for gauged U( 1) and Yukawa couplings, 
respectively, to a fermion ^ as in 

C C yip<p4>- (11) 


On the other hand, the scalar self-coupling also receives corrections like 

1 


Am —> Am + 


EY+z-E^ 1 ) 


( 12 ) 


v ' 647r 2 

with Y2 denoting the summation with respect to degrees of freedom. Therefore, in fact, the 
condition (J9J) can be satisfied only if accurate cancellations happen in Eq. f[T 2 ]) . Realization 
of a thermal inflation by the potential (J2]) is unstable against radiative correction, and 
in this sense, successful thermal inflation by this potential is difficult. This would be a 
reason why Lyth and Stewart have considered the flaton potential lifted by higher-order 
nonrenormalizable terms in the original “thermal inflation” paper 9J, as we will also do. 

Before turning to higher-order potential, we estimate the number of e-folds in the axion- 
dominated Universe by the dominated PQ scalar held decay, because the above Eqs. (ED, 
0, and © have been derived in a radiation-dominated Universe background. In the KSVZ 
model, the PQ scalar mostly decays mostly into two axions and a little into SM particles. 
We use Br to parametrize the decay branching ratio of the PQ scalar into SM particles. 
Then we have 


Ptotal PSMrad “T Paxiom 

PSMrad ■ Paxion Bl . 1 Bv, 
T 4 

PSMrad ' 


(13) 

(14) 

(15) 
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Since thermal inflation starts in the axion-dominated Universe, the condition ptotai — 
paxion — V{p = 0) is rewritten as 


1 

Br 30 ' 


1 /( 0 ). 


(16) 


We obtain 

T ‘ = \B? 30*) (W) 7 ’ (17) 

Ni = -n 4 ^-jln(^ 9 .)-jln(^:)-iln^. (18) 

In the axion-dominated Universe, the temperature interval between T t and Tf becomes 
shorter because the energy density of thermalized radiation is subdominant. As a result, the 
number of e-folds also becomes small, which can been seen as the effect of the third term 
in Eq. (TT8l) . We list various physical quantities about thermal inflation in Table H] One can 
easily find the hierarchy of inequality ([9]) for realizing even just N ~ 1. 


TABLE I: Quantities in thermal inflation by the potential d2]) 


A 

h 

v(GeV) 

Tj(GeV) 

7/(GeV) 

N 

A N eff 

Tr( GeV) 

10~ 6 

2.54 

10 8 

7.43 x 10 5 

2.73 x 10 5 

1.00 

0.05 

1.6 x 10 6 


B. Thermal inflation by a higher-power potential 


As we have seen, thermal inflation by a wine bottle potential can realize the expansion of 
only 0(1) e-folds if the condition ([9]) is satisfied. Hence, in this subsection, as a reference, we 
estimate the number of e-folds in the case where the potential is lifted by a nonrenormalizablc 
higher-order term. The scalar potential is expressed as 

I ( n I < ^ n 

V(p) = V 0 -m 2 M 2 +(19) 

Here, we stress that the absence of the |<^| 4 term is not necessary. We just assume that it is 
negligible due to the small coupling constant. Vo, the VEV, and the mass of cp are expressed 
as 


U 0 


n 


1) 


,2 n 


A 2 (n-2) ’ 


( 20 ) 
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( 21 ) 


, , /m 2 A 2 <"- 2 >\ 5 !^ 

”= <*>=(— 

y 2(n-l) 

— n(n — 1) ^ 2 (n-2 ) 5 (^ 2 ) 

by using the stationary condition and V(v) = 0. The number of e-folds by thermal inflation 
is give by 

N 2n = N 4 -~ In ?r + -(n - 2) In (—) , 

4 4(n — 1) 2 V ; \ v ) 

N 4 = — ln4v / 3- Inf—H—In - h, 

4 V 30 V 2 m p 

with Mp being the reduced Planck mass. We list various physical quantities about thermal 

inflation in Table HH Thermal inflation can be realized with more sensible coupling h than 

that in Table [B 


(23) 

(24) 


TABLE II: Quantities in thermal inflation by the potential (1191) 


A(GeV) 

h 

v(GeV) 

Ti( GeV) 

7) (GeV) 

N 

AN eff 

Tr( GeV) 

10 16 

8.27 x 10“ 3 

10 8 

2.79 x 10 3 

1.03 x 10 3 

1.00 

0.05 

5.9 x 10 3 

10 16 

8.27 x 10“ 2 

10 10 

2.79 x 10 6 

1.03 x 10 6 

1.00 

0.05 

5.9 x 10 6 

10 16 

8.27 x KT 1 

10 12 

2.87 x 10 9 

1.03 x 10 9 

1.00 

0.05 

5.9 x 10 9 


III. RELIC ABUNDANCES 

A. Axion dark radiation 

As we have seen, if the PQ scalar field dominates the energy density of the Universe, its 
decay produces many axions, and the Universe ends up with relativistic axion domination. 
When the total energy density p to tai from dominated axion p a and subdominant radiation 
Prad becomes comparable with V(tp), t = ti, the thermal inflation begins. Then we have 

= 0) = Ptotai (U) — "g^Prad(U) = ~j^Prad(tf) ^ a (t) ^ ' ( 23 ) 

Here tf stands for the time when the thermal inflation ends. After the thermal inflation, tp 
decays into SM particles and potentially non-SM particles again. Let us define the decay 
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branching ratio into non-SM particles br for later convenience. By using Eq. ([25]) . the ratio 
of axion to radiation is expressed as 


Prad | h =r = (1 — br)e LN Pa\H=V ( H r ^ , 


(26) 


where N and T are the number of e-folds of the thermal inflation and the decay rate of tp. 
The resultant axion dark radiation contribution is expressed in terms of AN e ff as 

43 /43/4\ 1/3 p a 


AN eff ~ ~ 


V .9* / 


x 


Prad 


H =r 


(27) 


B. Reheating temperature and possible baryogenesis scenarios 


The reheating temperature after inflation is estimated as 


T r = 


90 


-(1 - br)T 2 M\ 


1/4 


(28) 


N 2 g*{T R ) 

In the following, we estimate a reference value of the reheating temperature under the 
assumption of the instantaneous reheating T = H(tf ), which gives the highest reheating 
temperature. The availability of baryogenesis mechanisms depends on the reheating tem¬ 
perature after thermal inflation T r . 

For high enough reheating temperature Tr > 10 9 GeV, thermal leptogenesis by the 
lightest heavy RH neutrino decay of those with hierarchical masses is one of the simplest 
scenarios of baryogenesis [ll, 1^ ]. 

Nonthermal leptogenesis by RH neutrinos with hierarchical masses is available for a re¬ 
heating temperature 10 9 GeV > Tr > 10 6 GeV 13]. If this local 17(1) is in fact the gauged 
U(1)b~l symmetry, ip is identified with the Higgs held to break this symmetry with the 
B — L charge 2, and the decay <p into two RH neutrinos Nr is nothing but nonthermal 
production of Nr. 

For Tr < 10 6 GeV, low-scale thermal leptogenesis requires an enhancement of CP vi¬ 
olation. Here, for information, we refer to two examples. One possibility is the so-called 
resonant leptogenesis, where two RH neutrino masses are strongly degenerated and CP vi¬ 
olation is enlarged due to RH neutrino self-energy [yj]. Another way is an extension of the 
Higgs sector. It has been shown that, in the so-called neutrinophilic Higgs model, large 
enough CP violation can be obtained for the lightest RH neutrino mass of 0( 10 4 ) GeV 15]. 
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Another promising scenario would be electroweak baryogenesis. For a recent review, see, 


e.g., Ref. 




C. Results 


Here, by using Figs. |Tj and [2j we summarize the viable parameter space and available 
baryogenesis mechanism for some benchmark points we have studied. Inflation with the 
potential (( 2 ]) is labeled as n = 2, and n = 3 indicates inflation by the potential (flUl) with 
n = 3. Although we comment on cases with n — 2, one should remember that realization 
of thermal inflation by the n — 2 potential is, in general, difficult as stated above, so this 


is only for information purpose. In order to have a large enough CP violation e > 10 


-6 


m 


the Nr decay, which is equivalent to the so-called Davidson-Ibarra bound of the lightest RH 


neutrino mass for leptogenesis M Nr > 10 9 GeV 


3 , 0 , 


we suppose M Nr — 10 9 GeV and 


that the decay p —» NrNr is kinematically forbidden for m ^ < 10 9 GeV. We consider two 
cases of the PQ scalar VEV, v = 10 10 and 10 12 GeV. As far as implications for baryogenesis 
are concerned, the conclusion is the same for v < 10 10 GeV. 


1 . n = 2,v = 10 12 GeV case 

For most of the parameter space, we have Tr > 10 9 GeV. Thus, the standard thermal 
leptogenesis could work. 

2. n = 2 ,v = 10 10 GeV case 

For most of the parameter space, 10 9 GeV > Tr> 10 6 GeV is realized; however, m v < 10 9 
GeV. Nonthermal leptogenesis by the ip decay does not work as long as a hierarchical mass 
is assumed. Thus, a low-scale thermal leptogenesis with an enhanced CP violation or the 
electroweak baryogenesis with the extension of the Higgs sector is needed. 

3. n = 3 ,v = 10 12 GeV case 

For most of the parameter space, we have Tr > 10 9 GeV. Thus, the standard thermal 
leptogenesis could work. 
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4- n = 3, v = 10 10 GeV case 


As for the n = 2,v = 10 10 GeV case, Tr > 10 6 GeV is realized; however, rn v < 10 9 
GeV. Nonthermal leptogenesis by the ip decay does not work because the <p decay into them 
is kinematically forbidden. Thus, a low-scale thermal leptogenesis with an enhanced CP 
violation or the electroweak baryogenesis with the extension of the Higgs sector is needed. 


n=2,v=10 10 [GeV] 


n=2,v=10 12 [GeV] 




FIG. 1: Contours of the resultant AA e g = 1,0.4(thick red), 0.1, and 0.01 with solid lines from left 
to right, the mass of p> with dashed lines and the possible maximal reheating temperature after 
thermal inflation Tr with long dashed lines by the potential (f2|) . The shaded region corresponds 
to A N e ff > 0.4 which is disfavored by the Planck (2015) data. 


n=3,v=10 10 [GeV] n=3,v=10 12 [GeV] 




FIG. 2: The same as Fig. |T]but for the potential (1191) with n = 3. 
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IV. SUMMARY 


In this paper we have investigated scenarios with successful thermal inflation by a gauged 
1/(1) Higgs flaton held to dilute axions generated by late decay of the dominated PQ scalar 
held which dominantly decays into relativistic axions. We find the available parameter space 
to dilute axions. By estimating the flaton Higgs boson mass and the possible highest re¬ 
heating temperature after thermal inflation, if the U( 1) symmetry is the gauged U(1)b-l, 
we find that a promising viable baryogenesis in this cosmology is not nonthermal lcptoge- 
nesis with hierarchical RH neutrino masses but high- or low-scale thermal lcptogenesis or 
electroweak baryogenesis. 
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